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Abstract 

The nonexponential relaxation and aging inherent to complex dynamics manifested in a wide 
variety of dissipative systems is analyzed through a model of diffusion in phase space in the presence 
of a nonconservative force. The action of this force establishes a heat flow which maintains the 
system away from equilibrium. The inability of the system to find its equilibrium state becomes 
apparent through the presence of an effective temperature field. This is the temperature of the 
stationary nonequilibrium state reached by the system satisfying a generalyzed version of the 
fluctuation-dissipation theorem. The presence of a nonequilibrium temperature leads to a hierarchy 
of relaxation times responsible for the aging phenomena and to a relation similar to the Vogel- 
Fulcher-Tammann law. 
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I. INTRODUCTION 



Under certain conditions dissipative systems give rise to an often complex dynamics. This 
is the case of glassy systems which constitute an ubiquitous example of dissipative systems 
following a complex dynamics^. Glass is a nonequilibrium system whose dynamics is non- 
linear and has nonexponential relaxation functions which leads to the aging phenomena. 
Similar to glass is the process of protein folding which involves crossing a large number of 
energy barriers^i,— . Anomalous transport processes and strange kinetics based on fractional 
generalizations of the diffusion and Fokker-Planck equations^,— constitute other typical ex- 
amples. Hence, complex dynamics is the object of active research due to its multiplicity 
of manifestations. Here, we will focus on the nonexponential character of the relaxation 
functions and on aging in complex dynamics. The mechanisms underlying slow relaxation 
and aging still lack a clear and definitive elucidation. In glasses and proteins for example, 
several experiments and computer simulations have been done which support the explana- 
tion of these relaxation phenomena in the framework of the energy landscape paradigm as 
the result of activated diffusion through a rough energy landscape of valleys and peaks^. 

To contribute to the understanding of what these mechanisms are, we propose here a 
simple model to show a possible origin of nonexponential relaxation and aging based on the 
idea of the energy landscape and nontrivial energy barriers. This model which consists of 
the diffusion in phase space, provides a direct link between the phase space dynamics and 
the slow relaxation of the functions of the configuration of the system in the corresponding 
energy landscape. The slowing down of the dynamics clearly appears as a consequence of 
the freezing of some degrees of freedom which takes the system out of equilibrium. This fact 
is indicated by the presence of an effective temperature field incorporating the information 
of the suppressed degrees of freedom and depending on the equilibrium temperature at the 
moment the quench was applied. 

The paper is distributed as follows. In the second section we develop the thermodynamic 
framework in which our model is embedded. In section three, we write the Fokker-Planck 
equation governing the dynamics and study the relaxation. Section four is devoted to deriv- 
ing the generalized diffusion equation and obtaining the relaxation function and the hierarchy 
of relaxation times in some particular cases. Section five deals with the computation of the 
generalized fluctuation-dissipation theorem GFDT. Finally in section seven, we discuss our 
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main conclusions. 



II. THERMODYNAMIC FRAMEWORK 

We model the nonexponential relaxation phenomenon in a dissipative system as the Brow- 
nian motion of a particle of unit mass in a nonperiodic potential V{x), able to constitute a 
schematic representation of a rough energy landscape. Moreover, we will apply the methods 
of nonequilibrium thermodynamics at a mesoscopic level which consists of computing the 
entropy production from the Gibbs entropy postulate and postulating the linear phenomeno- 
logical equations relating the fluxes and the thermodynamic forces occurring therein. This 
constitutes a mesoscopic generalization of the macroscopic formalism of the nonequilibrium 
thermodynamics, the mesoscopic nonequilibrium thermodynamics (MNET )™ ,— . 

Thus, our thermodynamic study arises from the assumption that we are dealing with a 
physical ensemble of Brownian particles, a thermodynamic system to which thermodynamic 
relations can be applied^^.. We consider that the ensemble is initially at equilibrium in a 
heat bath at temperature To, and thus, is distributed according to the canonical distribution 



i^eq ~ exp 



,2 



(1) 



where H{T) = l/2u + V{x) is the Hamiltonian, and T = represents a point of the 

one-particle phase space, with x being the position and u the velocity. At a certain moment 
(t = 0), we apply a nonconservative force f{t) quenching the system. It is assumed that the 
nonequilibrium entropy of the system is given through the Gibbs entropy postulato^^,— 

S{t) = -k I ^\j,^dT + Seq, (2) 
J Weq 

where Seq is the equilibrium entropy. Variations in the probability density '?/'(r,t) imply 
changes in the entropy which can be obtained from Eq. These are given by 



(3) 



where we have taken into account that 6Seq = — J fieq/TodipdT, with figg = —V{x)+fio being 
the equilibrium chemical potential per particle of the Brownian gas and /Iq a constant . On 
the other hand, we assume local equilibrium in the phase space T consisting in postulating 
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the Gibbs equation at a local level in the phase spaced,—, which in the case of a homogeneous 
bath, reduces to 

T^8S = -j /i(r, t)5^{T, t)dT , (4) 

where the thermodynamically conjugated of the density ^/'(r, t) nonequilibrium chemical po- 
tential fi{T,t) appears. The expression of the nonequilibrium chemical potential is obtained 
by comparing Eqs. © and (jH), thus 

n{T,t) = kToln^- + fXeq . (5) 

A gradient of the chemical potential induces a diffusion process in phase space in order 
to restore the equilibrium state. Throughout this process, the distribution of probability 
changes according to the following Generalized Liouville equation 

-i;{r, t) + Vr{T, t) ■ Vr^(r, t) = -^^(r, t)f{t) - — j(r, t) , (6) 

where Vr(T,t) = {u, — is the velocity corresponding to the Hamiltonian flow, J(T,t) 
constitutes a diffusion current and Vr = V), where V stands for the spatial derivative. 
In order to completely characterize the diffusion process, one has to find the expression 
of the diffusion current by following the rules of nonequilibrium thermodynamicsi^. In 
this framework, through the computation of the entropy production one can establish phe- 
nomenological relations between currents and conjugated thermodynamic forces. Thus, from 
Eq. ^ we can determine the rate of change of the entropy 

To^ = -//i(r,t)^^(r,t)rfr, (7) 

which after substituting Eq. (jH)) and by using Eq. (jSJ becomes 

= {(fit) - vv(x)) - 1 j(r, t)|^/i(r, t)dY , (8) 

obtained after partial integration. Here, 

((/ W - V\/(x)) u)= I (fit) - \/V{x)) n^(r, t)dT (9) 



is the power supplied by the net force which is dissipated in the system as heat. Thus, the 
first term on the right hand side of Eq. (jHJ can be interpreted as the rate of heat exchange 
with the surroundings 

f^=-{{f{t)-VV{x))u) , (10) 



with dq being the amount of heat released in a time dt. However, the second term on the 
right hand side of Eq. (jH)) constitutes the entropy production a > 0, according to the second 
law, containing the current and its conjugated thermodynamic force {d/du)fi{r, t), which by 
postulate of the nonequilibrium thermodynamics are related through the phenomenological 
law 

J(r,t) = -^|-/i(r,t), (11) 

where L is a phenomenological coefficient. Therefore, Eq. (jHJ may be rewritten in a more 
significant form 

dS 1 dq , , 

expressing the entropy balance between the exchange of heat with the surroundings and the 
entropy generated in the irreversible processes established in the system when it is removed 
from its equilibrium state. By using the expression of the nonequilibrium chemical potential 
(0) in (fTTj) . one obtains 

J{T,t) = -j(^kBTo-^ + u^ij , (13) 

where we have performed the identification L/ipTo = 7, with 7 being the friction coefficient 
of the Brownian particle so that 



a 



7 f{kBToid/du)^{r,t) + ^{r,t)uf ^^ ^^^^ 



To J ^{T,t) 

One can conclude that is possible the existence of a stationary nonequilibrium state {dS/ dt 
0) in the system, i.e. not given by an equilibrium maxwellian, for which 

1 dstq 



Tn dt 



-(Tst . (15) 



III. FOKKER-PLANCK EQUATION AND RELAXATION 

The Fokker-Planck equation governs the dissipative dynamics when the description of 
the state of the system is given in terms of the probability density ip{T, t). This equation is 
obtained after substituting the expression of the current (fT^ into Eq. (jH)) leading to 

g^V'(r,«) = -g^wP(T,t) + g^iAff.i) (VV-(rc) - /(«)) (16) 
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We will assume that u is the fast variable, thus the dynamical processes in the system are 
associated to configurational changes related to x, the slow variable. Thus, we write 



^/'(r,t) = 0,(u,t)p(a;,t) , (17) 

where t) is the conditional probability . The probability density is p{x, t) = J 'ipiT, t)du 
evolves according to 

— p(x,t) = -V / uij(T,t)du , (18) 

obtained by partial integration of Eq. p6|) . which defines the current J{x, t) = J wp{r, t)du. 
This current satisfies the equation 

^J(a;,t)+7J(x,t) = (19) 
- {p(x, t) {VV{x) - fit)) + ksVpix, t)T{x, t)} . 

which can be obtained after multiplying Eq. (fT^ by the velocity u and performing the 
partial integration, by using the decoupling (fTTj). Here, kBT{x,t) = f u^(j)x{u,t)du is the 
second moment of the conditional distribution (f)x{u,t)^ playing the role of an effective tem- 
perature which contains information on the frozen degrees of freedom. This temperature 
field introduces thermal barriers in the systen>^. The presence of these thermal barriers is 
a nonequilibrium effect that disappears when the system is at equilibrium. Over long time, 
when the transient regime has died out, Eq. (|T9|l becomes 

J(x, t) = -ro{p(x, t) {VV{x) - fit)) (20) 
+ kBVp{x,t)T{x,t)} 

where tq = 7"-^. 

On the other hand, in Eqs. (fT^ and the current J(x, t) appears coupled to the 
temperature field T{x,t) which evolves according to 

^p(x, t)T{x, t) = -^Vp(a;, t)h{x, t) (21) 

-2 (VV{x) - fit)) - -^p(x, t) (T(x, t) - n) , 

Kb KbTq 

obtained by multiplying Eq. by and integrating in u. Here, we have defined the heat 
flow 

h{x,t) = ^J u^(l),{u,t)du . (22) 



For times t S> tq, Eq. (PT|) reduces to 

p{x, t)kB iT{x, t) - To) = -ToVpix, t)h{x, t) (23) 
-To{VV{x)-f{t))J{x,t) . 

In the particular case of equilibrium, (j)x{u,t) would be an equilibrium Maxwellian, thus 
J{x,t) = h{x,t) = which would lead to T{x,t) = To. If, on the other hand, we integrate 
Eq. we obtain 

fc^(T(t)-T,) = ro^, (24) 

with T{t) = {T{x,t)). This last result together with Eq. (|Tn|l leads to the interesting 
conclussion that even if f{t) = the partial elimination of the fast variable would introduce 
a nonequilibrium temperature!^. Only for the particular case of a soft energy landscape with 
small gradients W{x), the difference T{t) — To would be negligible in comparison to T^. We 
can make this statement more precise by defining a critical length scale L = AV/W{x), 
where AV" is the characteristic variation of the potential, such that if {'jL)'^ /ksTo <^ 1, then 
T{x,t) = To. 

In fact the Brownian motion as described by the Klein-Kramers equation ()lfi|l is a prob- 
lem with multiple time scales as was pointed out in refs.-^,-^ and implicitely included in 
the generalization of the Smoluchowski equation found previously in^SJ^. Here, we have 
embedded all of these scales in the nonequilibrium temperature T{x,t) or T(t), constituting 
a kind of mean field approach which leads to the nonexponential relaxation and aging we 
will show in the following sections. 

To sum up, in this section we have obtained the evolution equations for the first moments 
of the conditional probability density which are of interest for the thermodynamic description 
of the Brownian gas. More detail and some consequences will be shown in the following 
sections. 

IV. GENERALIZED DIFFUSION EQUATION 

Here, we will establish the generalized diffusion equation and examine its consequences 
in respect to the relaxation phenomena. 
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After defining the effective potential t) = V{x)+kBT{x, t), Eq. (pUj) can be expressed 

J(x, t) = -D{x, t)Vp{x, t) - rop(x, t) (V$(x, t) - f{t)) (25) 
= Q{x,t) + Top{x,t)f{t) , 

where D{x,t) = TokBT{x,t) is the generahzed diffusion coefficient and Q{x,t) corresponds 
to the current of the purely relaxational system {i.e. in absence of nonconservative forces). 
Thus, by substituting Eq. into Eq. (fTBj) and assuming tq = 1 {i.e. rescaling the time 
t), this last equation becomes the generalized diffusion equation 

^p(a;, t) + V p(x, t)f{t) = VQ{x, t) . (26) 

In addition, as a consequence of the elimination of degrees of freedom, the dynamics becomes 
non-Markovian and depends on the equilibrium temperature Tq at the quench time. In the 
frame of reference defined through the transformation x* = x — J* f{v)dv, the current ()25|) 
reduces to Q{x*,t). Thus, in this frame one has a purely relaxational process for which 

^^p{x\t) = VQ{x\t) . (27) 

This equation admits a quasi-equilibrium solution 

Pqe (^^5 ^) 

exp | — /'^ ksTix' for which Q{x,t) = 0, where to simplify the nota- 

tion we have omitted the superscript * on the x. 



A. Population dynamics and relaxation function 

To examine the problem of the population dynamics with its origin in the generalized 
diffusion equation (|27p and to find the corresponding relaxation function we will assume that 
at a point X2 there is a sink where p{x2, t) = and a source at xi. Hence, a quasi-stationary 
current, Q{t) can be established in the system, thus, integration of Q{x,t) gives us the 
quasi- stationary probability density in terms of Q{t) 



On the other hand, a second integration leads to 

rx2 

D{y,t)pi,ge.{y,t) 



n{t) = Q{t) I dxp,e{x,t) r^, , (29) 



with n{t) = J^^ p{x, t)dx being the population at the left of the sink. Thus, after integration 
of Eq. (fTH)) we obtain 

jnit) = - {g(x,, t) - g(xi, t)} = -Q{t) , (30) 

which leads to the following rate equation 

d 



^^n{t) = -K{t;Xb)n{t) , (31) 



where 



K{t;xb)'^= / dxpge{x,t) / — - (32) 

A'l Jx D{y,t;to)pi.qe.{y,t) 

is the rate constant, and Xb is the position of the barrier between xi and X2- The relaxation 
equation (p?T| admits the solution 

n{t) = n{to) exp | - ^ K{t'; Xb)dt'^ (33) 

= n{to) exp{-g{t)} 

or 

n{t) = n{to) exp{-g{t/ro)} , (34) 

where n(t) = n(t/To). Both Eqs. (jH^ and (jMj) constitute an interesting result since we 
obtain a hierarchy of relaxation times starting from a Markovian equation, the Fokker-Planck 
equation (fTHjl . Here, the function g{t/To) might be interpreted as an algebraic function 
{t/r^Y obtaining a stretched exponential behavior^^ or as a logarithmic function —a log At 
which leads to a power law behavior^*^ that characterizes anomalous diffusion^^. In addition, 
another reading of our result is possible because assuming the form of g{t/TQ) is equivalent 
to assuming the form of the distribution of residence times il){t) = —-^n(ty^. Thus, an 
in-depth analysis of the implications of the supposition of a nonexponential distribution 
of residence times usually performed in Continuous Time Random Walk models, reveals 
that this hypothesis might also be rooted in the nonequilibrium character of the dynamics 
inherent to the energy landscape picture as we have shown. 

B. Hierarchy of relaxation times 

Inherent to the nonexponential relaxation Eqs. ()33|) and ()34|) there is a hierarchy of 
relaxation times. To show this clearly here, we will derive the hierarchy of relaxation times 
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for the particular case of a bistable potential and the state of quasi-equilibrium. Thus, let 
us assume here that V{x) is a bistable potential with one well much deeper than the other. 
In this case, we can obtain the dynamic of the population of the shalower well through an 
equation analogous to Eq. (j3H) . where now according to the Kramers formula 



(35) 



where we have defined the frequency /o = ^■\/V"{xi) \ V"{xb)\, and the barrier height A = 
V{xb) — V{xi). Thus, 

n{t) = n{to)exp{- / > (36) 



I. Jto 

with 

From Eq. ()24j] . we can write 

T(t)-T(to)=/(tf -/(^o)^ (38) 
where now T{t) = {T{x,t))^^ and (VV {x)u) = 0. Therefore, 

or equivalently 

where T*{t) = T{t) — To, i.e. a scale of temperature whose zero point coincides with the 
temperature of the heat bath. Here, assuming that — T{to) ^ 1 we can neglect the 
second exponential on the right hand side of Eq. pn|) , leading to the result 

similar to the celebrated Vogel-Tammann-Fulcher equation. In the next order of approxi- 
mation we would have 

providing an expression for the relaxation time which is finite when T*{t) = T*{to). 
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V. FLUCTUATION-DISSIPATION THEOREM 



Here, we want to derive the fluctuation-dissipation theorem for the nonequihbrium sys- 
tem, which is locally satisfied and provides us with a way to compute the nonequihbrium 
temperature. In order to highlight the main details and avoid lengthy calculations, we 
will homogenize Eq. ()27|) introducing a preaveraging condition with respect to the quasi- 
equilibrium distribution 

l'' = |:{^<"^''<^-" + ''<^-*'l*<^--')} • 

where D{t) = kBT{t) is the preaveraged diffusion coefficient. Thus, now the 
quasi-equilibrium density becomes a local Boltzmann probability density pqe{x,t) ~ 
exp{-^x,t)/kBT{t)} . 

Consider a time- dependent external field F{t) = FQ{t — t^,) applied to the system Eq. 
()43|) at the waiting time when this system is in the quasi-equilibrium state described 
by pqe- Here, Q{t — tw) is the step function. This force modifies the potential, and we 
assume that the perturbed potential is t) = $(a;,t) — F{t)B{x). If the field is weak, 
the deviation with respect to stationarity of the average value of any function A{x) of the 
configuration of the system is a linear functional of the field 

{A{t))^ - {A{t))^^ = f R{t, t')F{t')dt' , (44) 

which introduces the response function R(t,t'). Now {A(t))p can be calculated once the 
perturbed distribution function t) is known: 

{A{t))p = J A{x)p{x,t)dx . (45) 

Since the system is in the nonequihbrium stationary state at the time t^, p{x,t) is related 
to Pqe{x,tw) thrugh the perturbed Green function G{x,t \ x',t^) 

p{x,t) = j G{x,t I x',t^)pqe{x',t^)dx' , (46) 

and 

exp(Fi?(x)A^T(t)) 
^^^^"^'^^ {exp{FB{x)/ksTm/'^^''''> ' ^^^^ 

To compute {A{t))p by using Eq. ()45|) with the help of Eq. ()46p we have to know the relation 

between G{t,x \ x',tuj) and the corresponding unperturbed Green function G{t,x \ x',tw)- 
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This relation can be established from the identity 

Pqe{x,t) = / G{t,X \ X',ty,)pge{x',t^)dx' , (48) 



by applying the relation (jT7|) at both sides of Eq. (PS|). leading to 

CIt X I X- t ) exp(Fi;(x)/t,r(t„,)) g {e.p{FB(x)/ksTm„ 
^ ■ ' ' " exp (FB{x)/kBm) ''I'' ^ I ^ ■ (exp (FB(x)/feT((„))),, ' 

By using Eqs. (US)), (jlH)) and (UHl), after expanding the exponentials with respect to F, one 
then obtains 

(Am), - = F {^C,«(M) - ^C,.«(Mj} . (50) 

where CA,B{t,t') = {A{t)B{t')) - {A{t)) {B{t')). We assume that the state of the system 
varies slowly enough so that it is safe to replace T(t) for T{t^). Thus, one achieves 

m)), - {A{t))^^ = ^-i^ |7CA,B(t, t')F{t')dt' , (51) 

which defines the nonequilibrium response function 

Rit, t^) = -l—^CA,B{t^ tu,) , (52) 

for t > tw, constituting a local version of the fluctuation-dissipation theorem, the GFDT. 

This section shows that the kind of effects related to the nonequilibrium nature of the 
solutions araised in the Fokker-Planck equation can be tested by means of linear response 
theory. 



VI. CONCLUSIONS 



Here, we have learned that the dynamics of a brownian gas under the action of a noncon- 
servative force slowly varying over time is representative of the complex dynamics observed 
in dissipative systems. After quenching, the system is taken out of its initial equilibrium and 
then a heat flow is established. We model the effects of the quench of the system through 
the elimination of the fast degrees of freedom from the dynamics described by the Fokker- 
Planck equation. The departure from equilibrium and consequently the exchange of heat is 
maintained by means of the nonconservative force. We have also considered the existence of 
a nonperiodic ruogh potential in order to incorporate activated processes into our analysis 
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according to the energy landscape paradigm inherent to several complex systems. A series 
of facts constituting the main results of our analysis follows. 

First, the frozen degrees of freedom are manifested in the constrained dynamics through 
an effective temperature field different from the equilibrium temperature. The average of this 
temperature field plays the role of the temperature that characterizes the non-Markovian 
diffusion equation obtained from the Fokker-Planck equation, which describes the dynamic 
of the slow degrees of freedom . We also derive the relaxation equation for the effective 
temperature field from the Fokker-Planck equation which enables us to give a sound base 
to the nonequilibrium temperature as a result of the balance between the added heat and 
the heat generated during the relaxation process in the system. In addition, we show that 
the nonequilibrium temperature computed using the quasi-equilibrium distribution, i. e. the 
stationary solution of the non-Markovian diffusion equation, characterizes the generalized 
fluctuation-dissipation theorem we have derived, thus coinciding with the effective temper- 
ature referred to by several authors in the literature^,—. 

On the other hand, we find a hierarchy of relaxation times as a consequence of the 
non-Markovian character of the constrained dynamics of the system. In this respect, our 
model also describes stretched exponential relaxation as well as anomalous diffusio n^^,^^ 
and other Continuous Time Random Walk problems^^. Likewise, as a consequence of the 
nonexponential relaxation, we obtain an equation similar to the celebrated Vogel-Tammann- 
Fulcher law which is, in our opinion, an important result. 

Other approaches based on the energy landscape, considering the relaxation as a series 
of activated processes, assume that the activation energies are independent from the tem- 
perature but with a certain random spatial distribution. Averaging over this distribution of 
energy barrier gives one the corresponding relaxation^,—,—. In the literature another types 
of models undergoing nonexponential relaxation can be found, as for example, disorder 
modela^^. 

We conclude that our model incorporating the multiplicity of time scales embedded in 
the Fokker-Planck equation, although simple, contains the main features in the dynamical 
slowing down observed in a wide variety of complex systems. On the other hand, with 
our approach we reach more insight avout the physical meaning of the corrections to the 
Smoluchowski equation found in previous works as those referred to in this contribution, 
stating this problem in the modern language of the complex dynamics. In addition, our 
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model constitutes a way to study Brownian motion in disordered space leading to results 
similar to those reached through the approach based on fractal Fokker-Planck equations 
. Therefore, we think this work will contribute to the understanding of some aspects of 
complex dynamics in several dissipative systems. 
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